Q-ball formation in the gravity-mediated SUSY breaking scenario 
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We study the formation of Q-balls which are made of flat directions that appear in the supersym- 
metric extension of the standard model in the context of gravity-mediated supersymmetry breaking. 
The full non-linear calculations for the dynamics of the complex scalar field are made. Since the 
scalar potential in this model is flatter than <j!>^, we have found that fluctuations develop and go 
non-linear to form non-topological solitons, Q-balls. The size of a Q-ball is determined by the most 
amplified mode, which is completely determined by the model parameters. On the other hand, 
the charge of Q-balls depends linearly on the initial charge density of the AfHeck-Dine (AD) field. 
Almost all the charges are absorbed into Q-balls, and only a tiny fraction of the charges is carried 
by a relic AD field. It may lead to some constraints on the baryogenesis and/or parameters in the 
particle theory. The peculiarity of gravity-mediation is the moving Q-balls. This results in collisions 
between Q-balls. It may increase the charge of Q-balls, and change its fate. 
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I. INTRODUCTION 

A Q-ball is a kind of a non-topological soliton, whose 
stability is guaranteed by some conserved charge in scalar 
field theory [0j2j. It can be made of the scalar fields 
which appear as fiat directions in the supersymmetric ex- 
tension of the standard model |^,^. Particularly, in the 
minimal supersymmetric standard model (MSSM), the 
baryon and/or lepton number are the conserved charges, 
since those fiat directions consist of squarks and/or slep- 
tons It is known that large Q-ball solutions exist 
when both gauge-mediated and gravity-mediated super- 
symmetry (SUSY) breaking scenarios are included 
In the gauge-mediation scenario, the baryonic charged 
Q-ball, the B-ball, is stable against decay into nucleons, 
since the energy per unit charge becomes less than the 
nucleon mass, 1 GeV, for large enough Q-ball charge: 
E ~ mQ3/4 g Therefore, large B-balls can be a 
promising candidate for the cold dark matter. On the 
other hand, Q-ball energy grows linearly in the gravity- 
mediation scenario: E mQ [^j. They can thus decay 
into both nucleons (baryons) and lightest supersymmet- 
ric particles (LSPs), which become the dark matter in the 
universe. In the both scenarios, we can expect a close re- 
lation between the energy density of the baryon and dark 
matter such as J7b ^ ^dm 1^,11 i^b and floM are density 
parameters of the baryon and the dark matter, respec- 
tively). In particular, a somewhat more definite relation 
on the number densities hold for the gravity-mediation 
scenario: ulsp — Nsfsnb |||,^- Here Nb is the number 
of LSP decay products from the scalar field (fiat direc- 
tion) with unit baryon number, and fs is the fraction 
of the charge stored in the form of Q-balls. For these 
mechanism to work, the charge of B-ball should be in 
the range lO^o - lO^o 

Those large Q-balls are expected to be created through 



Affleck-Dine (AD) mechanism in the inflationary uni- 
verse 1^-^. The coherent state of the AD scalar field 
which consists of some flat direction in MSSM becomes 
unstable and instabilities develop. These fluctuations 
grow large, and are expected to form into Q-balls. The 
formation of large Q-balls has been studied only linear 
theory analytically [^-^ and numerical simulations was 
done in one-dimensional lattices [1). Both of them are 
based on the assumption that the Q-ball configuration 
is spherical so that we cannot really tell that the Q-ball 
configuration is actually accomplished. Some aspects of 
the dynamics of AD scalar and the evolution of the Q- 
ball were studied in Ref. but the whole dynamical 
process was not investigated, which is important for the 
investigation of the Q-ball formation. 

Actual Q-ball formation is confirmed in our recent 
work where we showed the formation of Q-balls in 
the gauge-mediated SUSY breaking scenario using lattice 
simulations in one, two, and three dimensions in space. 
In that scenario, the typical size of Q-balls is determined 
by that of the most developed mode of linearized fiuctua- 
tions when the amplitude of fiuctuations grow as large as 
that of the homogeneous mode: {S4>^) ^ (j)^- Almost all 
the initial charges which the AD condensate carries are 
absorbed into the formed Q-balls, leaving only a small 
fraction in the form of coherently oscillating AD conden- 
sate. Moreover, the actual sizes and the charges stored 
within Q-balls depend on the initial charge densities of 
the AD field. We also find that the evolution of the Q-ball 
crucially depends on its spatial dimensions, and the sta- 
ble Q-ball can exist only in the form of three-dimensional 
object. 

One may wonder if these results are peculiar to the 
gauge-mediation scenario which has a very fiat scalar po- 
tential for the large field value. For a very fiat scalar 
potential, larger Q-balls are easily formed, because the 
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energy of the Q-ball grows E mQ^/^ Q: the larger 
the charge is, the smaller the energy per unit charge is. 
On the other hand, the Q-ball energy grows linearly in 
the gravity-mediation scenario: E ^ mQ j^. Thus, we 
naively expect less effective Q-ball formation, particu- 
larly for large charge Q-balls to form. 

In this paper, we show the Q-ball formation in the 
gravity-mediation scenario by the use of numerical cal- 
culations. We find it very similar to gauge-mediation 
version, but some different new features are revealed. 

In the next section, we see the origin of the fluctua- 
tions of the complex scalar field, and show the instability 
ban d. R esults from numerical calculations are shown in 
Sec. Here the charge and the size of Q-balls are found. 
In Sec. we will make some comments on the B-ball 
baryogenesis. We will show some peculiar phenomena of 
the Q-ball in the gravity-mediation scenario, such as the 
moving Q-balls, and their collisions as a result. Section 
VI is devoted to our summary and conclusions. 



II. INSTABILITIES OF AFFLECK-DINE 
CONDENSATE 

Q-balls with large charge are expected to be formed 
through Affleck-Dine mechanism It is usually con- 
sidered that the AD fleld are rotating homogeneously 
in its effective potential to make the baryon numbers. 
However, if we consider the SUSY-breaking included po- 
tentials, spatial instabilities of the AD field are induced 
by the negative pressure because of its potential being 
flatter than 0^ @JlillL To be concrete, let us take the 
following potential | : 
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where $ is a complex scalar field which brings a unit 
baryon number, A is a coupling constant of order unity, 
H is the Hubble parameter, c is a positive order one 
constant, M is a large mass scale which we take it as 
~ 2.4 X 10^^ GeV, and the ii'-term is the one-loop cor- 
rections due especially to gauginos, and the value of K 
is estimated in the range -0.01 to -0.1 In this 

potential, the pressure is estimated as |^] 
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where is the energy density of the scalar field (Here 
we assume that \K\ ^ 1 so that the first term in Eq.P) 
can be approximately rewritten in the power-law 0^+'^^. 
Therefore, the negative value of K is the crucial point for 
instabilities. 

The homogeneous part of the field evolves as 
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where we define the field $ to be 



(4) 



Then the equations for the linearized fluctuations can be 
written as 



3H6(I) ■ 



^V^6(l) - 29{t)(l){t)S9 ~ 9^{t)6(l) 



= 0, 



(l){t)6e + 3H[9{t)6(j) + (f>{t)S9] 



V^S9 



+2^{t)59 + 29{t)5^ = 0, (5) 

We are now going to see the most amplified mode. To 
this end, we take the solutions in the form 
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If a is real and positive, these fiuctuations grow exponen- 
tially, and go non-linear to form Q-balls. Putting these 
forms into Eqs.(||), we get the following condition for the 
non-trivial 5(j)Q and 59^, 



FiH) +0 + 0^ + — + 3m^K 
a-' 
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d + (i^ H ^ 
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= 0, 
(7) 

where F(H) = - - ^H^. 

It is appropriate to assume that H <^ m and a <C 
a, since the AD field oscillates when H ^ in, and the 
adiabatic production of fluctuations will occur. Then, 
Eq.(0) will be simplified as 
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Since 9^ ~ m^, for a to be real and positive, we must 
have 



3m^K < 



(9) 



As we are considering to be a negative value, an in- 
stability band will exist. This is because the oscillating 
scalar field in the potential flatter than has negative 
pressure, which leads to the instability of the homoge- 
neous field. Thus, the instability band should be in the 
range 



< ^ < 3m^\K\ 



(10) 
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We can easily derive that the most amplified mode is 
the center of the band: [kmax/aY — im'^\K\/2, and it 
corresponds exactly to the Q-ball size which is estimated 
analytically using the Gaussian profile of the Q-ball 
We will see shortly that it also coincides with the size 
actually observed on the lattices in our simulations. 



III. CHARGE AND SIZE OF Q-BALLS 

In this section, we show the results of the lattice sim- 
ulations. In the potential the AD field obeys the 
equation 



l + K + K\og 



o\2 

+ — |<i>|4<i> = 0. (11) 

Here we have calculated in the matter-dominated uni- 
verse, so that H = 2/'it. In the context of AD mecha- 
nism for baryogenesis, the A-terms, such as VA-term ^ 
{AX/M)(j)'^ + h.c., should be added to the potential (|l|) in 
order to make the AD field rotate around in its potential. 
Instead, we take ad hoc initial conditions and neglect A- 
terms, since they do not affect the later dynamics of the 
field crucially. Therefore, the AD field possesses some 
initial charge density. 

It is more convenient for numerical calculations to take 
the real and imaginary decomposition $ = (^i +i(f>2)/V2 
and rescale as follows: 
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For the initial conditions, we take some large vev in the 
real axis and put some angular velocity to the imaginary 
part. In addition, we put initial fluctuations very small 
values O(10~^). Thus, they have the form 



^i(o) = A + M(0, ^'M = m{0, 

^2(0) = 5C(O, if'M = D + 5D{0, 



(13) 



where A and D are some constants, independent of 
the position is space, 5 A, SB, 6C, and 6D are ^ depen- 
dent small random variables, and the prime denotes the 
derivative with respect to r. Notice that the important 
features of the dynamics of the field are not affected by 
how we take these random variables, if we do not choose 
very peculiar distributions. 

We have calculated the dynamics of the AD scalar field 
for various parameters, and find that the initially (ap- 
proximately) homogeneous AD field deforms into a lot of 
clumpy objects. See Figs. |l| and ^j. All of them conserve 
their charge very well, so they must be Q-balls. (We 
observed charge loss and exchange between two Q-balls 
in some cases. We will discuss them in Sect. 0.) The 
profile of the Q-ball is a spherically symmetric thick-wall 
type, and fits very well to the Gaussian. In these figures, 




FIG. 1. Configuration of Q-balls on three-dimensional 
lattice. More than 40 Q-balls are formed, and the largest one 
has the charge with Q ~ 5.16 x 10^^ 




FIG. 2. Configuration of Q-balls on three-dimensional 
lattice. In each direction, the box size is half of that in Fig. ^. 
More than 10 Q-balls are formed, and the largest one has the 
charge with Q ~ 1.74 x 10^*^ 



we take (pi{Q) — ¥'2(0) — 2.5 x 10^ for the initial con- 
ditions on the 64'^ three-dimensional (3D) lattices with 
= 0.1 and = 0.05 for the large and small lattice 
boxes, respectively. It seems that there is no box-size 
effects in these parameters, since these two figures look 
the same. They have similar charge distributions and the 
Q-ball size is the same, as expected from the analytical 
estimate, Rphys ~ \K\~^/^m~^ . Actually, the numbers 
of Q-balls with the charge larger than 10^^ are 7 and 2 
in the large and small box, respectively. 

Comparing to those Q-balls which appear in the gauge- 
mediated SUSY breaking scenario, the size of the Q-ball 
is much smaller for the same charge, and most of the Q- 
balls has the same order of size. This is because Rphys ~ 
|X|-i/2„|^-i fQj. ^Y^Q gravity- mediation, which does not 
depend on the charge Q, while Rphys m~^Q^/^ for the 
gauge-mediation. We thus observe large-charged Q-balls 
with relatively small size. 
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FIG. 3. Amplitude of the AD field after formation of 
Q-balls. This configuration is the slice at z — 6.3. The am- 
plitude of relic field outside the Q-balls is two or three orders 
smaller than that of the center of the Q-balls. 



As in the case of the gauge- mediation scenario , we 
observe almost all the charge which initially AD conden- 
sate has absorbed into Q-balls, and the amplitude of the 
relic AD field is highly damped. This means that the 
fraction of the charge outside Q-balls is very small. Fig- 
ure ^ shows the amplitude of the AD field of the slice at 
z = 6.3 in the larger box for another realization of sim- 
ulations. Notice that there is relic field outside Q-balls, 
but the fluctuations are rather large, and we may not 
be able to consider it as a homogeneous condensate. In 
particular case of Figs. |l| and || Q-balls carry more than 
97% and 99% of the total charge, respectively. In Fig. ^, 
the fraction of the charge outside the Q-balls is shown as 
a function of the number of Q-balls which we take into 
account. In the larger box simulation, only seven of the 
largest Q-balls hold more than 95% of the total charge. 
On the other hand, more than 97% is stored in only two 
of the largest Q-balls in the small box one. Notice that 
the dotted line (small box) is below the solid line (large 
box), because the resolution is twice as good in the for- 
mer simulation: the lower bound is determined by the 
resolution of each simulations. 

Analytically, some features of the Q-ball in gravity- 
mediation is known M. For example. 



E ^ mQ, , R 
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-1/2, 



m, etc. (14) 



They are all confirmed numerically. One example is 
shown in Fig. |^. This confirms the first relation of 
Eg. (p^ , which implies that the energy per unit charge 
is constant of 0{m). 

It is the best way to investigate the dynamics of Q-ball 
formation on f/iree-dimensional lattices, but it is practi- 
cally difficult to do, since we need somewhat high res- 
olution, and many runs for various parameters to look 
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FIG. 4. Fraction of the charge outside the Q-balls. The 
solid and dotted lines denote the results from the simulations 
shown in Figs. |l] and respectively. 
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FIG. 5. Dependence of the energy of the Q-ball on its 
charge calculated on three-dimensional lattices. This confirms 
the analytical estimate: E ~ mQ (the dotted line). 



at. Thus, we also calculate on one and two-dimensional 
lattices for more rigorous quantitative analysis. There- 
fore, we must know the evolution of Q-balls after their 
formation. We follow the similar discussion we made for 
Q-balls in the gauge- mediation scenario ||l^]. Since a Q- 
ball configuration is the energy minimum with some fixed 
charge Q, Q is constant with respect to time, so 



Q = a^Qo ~ a^R^q ^ const. 



(15) 



where Qd the charge in D dimension, and q = 4>i(f>2 — 
(j}i4>2 is the charge density. If we assume the form of a 
Q-ball as 



(/)(x, t) = 0(x) exp(iwt), 
the energy of a Q-ball can be calculated as 



(16) 
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where 



grad 
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(18) 



Here we assume that the logarithmic term of the first 
term in the potential (|l|) is small compared to the unity, 
so that we can approximate it in the power-law form. 

When the energy take the minimum value, the equipar- 
tition is achieved: Egrad ~ Vi and Eg^ad ~ V2. From 
these equations addition to the charge conservation, we 
obtain the following evolutions: 



JZ-D)\K\/[1+{D-1)\K\] 



(19) 



which we observed approximately the same features nu- 
merically. For = 3, we get very natural relations: 
Rphys = Ro. ~ const., Lo ^ const., and <j) ^ const. Al- 
though (j) decreases as time goes on for D = 1 and 2, R 
and u! is almost constant, since \K\ <C 1. This feature is 
different from that in the gauge mediation scenario, and 
is good for long simulations because low-dimensional Q- 
balls do not shrink the size so much. 

Now we will see that the size of the Q-ball is deter- 
mined by the most amplified mode. Comparing to the 
actual sizes observed on lattices, we also calculated nu- 
merically for linearized fluctuations. Although we decom- 
posed the complex field in radial and phase direction in 
the previous section, it is more convenient to decompose 
it into real and imaginary part for numerical simulations. 
We thus integrated the following mode equations in di- 
mensionless variables: 



Sip'- + 3hip'^ + 



fc2 
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Ifl+Lpl 
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Sipi 

6ipj = 0, (20) 



where {i,j) — (1,2), (2,1), and to = m/M. 

Figure ^ shows the power spectrum calculated from 
a lattice simulation and the above linearized equations 
at T = 5.5 X 10^ and r = 6 x 10^. We take the lat- 
tices with lattice size N = 1024 and lattice spacing 
= 0.1 in one dimension here, because we need high 
resolution data to make the power spectrum smooth for 
lower k. These two times are just before and after the 
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FIG. 6. Power spectra of fluctuations of AD scalar fleld 
{k^\S(pk\'^ , D — 1) when the amplitude of fluctuations be- 
comes as large as that of the homogeneous mode: {Sip ) ^ ip . 
The top panels (a) and (c) show the full fluctuations calcu- 
lated on one dimensional lattices, while the bottom panels 
(b) and (d) show the linearized fluctuations without mode 
mixing. 



fluctuations are fully developed: {SpP') ^ pP' . For lin- 
earized fiuctuations, the instability band is almost the 
same as Eq.(^o|). For example, the upper bound is esti- 
mated by fc/TO = \/3a(T)|i^|i/2 ~ 2.5 for \K\ = 0.01 and 
T — 5.5 X lO'^. See panel (b). Even before the full de- 
velopment of fluctuations (panel (a)), rescattering effects 
kick the lower mode to higher, and the spectrum gets a 
little broader |l^ . Needless to say, the spectrum becomes 
extremely broad and smooth after fluctuations are fully 
developed (panel (c)). At any times, however, the peeks 
are at the same points for both lattices and linearized 
cases, and correspond to the typical size of Q-balls actu- 
ally observed on the lattices. Therefore, we can conclude 
that the size of the Q-ball is determined by the most am- 
plified mode of the linearized fluctuations when they are 
fully developed. For the case of Fig. ||, the typical size is 
^max ^ 0.5, which implies Rphys ~ a{Tf)/kmax ^ 28.9, 
where t/ ~ 5.5 x 10^ is the formation time. This value 
exactly coincides with the sizes of Q-balls observed on 
three-dimensional lattices. Actually, they are (a few) x 10 
in the dimensionless units. 

The actual values of the charge depend on the val- 
ues of the charge density which AD field initially pos- 
sesses. Since initial charge density is written as g(0) — 
(/Ji(0)(/92(0) for our initial conditions, we must check the 
dependence on both initial amplitude ipi (0) and angular 
velocity </32(0) of AD field. Results are shown in Fig. |^. 
Here we plot the largest charge Qmax against the initial 
AD charge density q{Q). We investigate two situations. 
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FIG. 7. Dependence of charges on the initial charge 
density g(0) = (^1(0)932(0) carried by the AD condensate 
on one-dimensional lattices. Open squares denote the case 
(/3i(0) = ¥'2(0), "pluses," crosses, and solid triangles denote 
the dependence on ^2(0) with 1/^1(0) fixed at 10^, 10®, and 
10^, respectively. 



The first one is changing both equally while fixing the re- 
lation </'i(0) = </'2(0), which is shown by open squares in 
the figure. This corresponds to the "maximum charged" 
Q-balls in terms of Ref. . We can fit all of these on 
the straight line (dotted line), Qmax ~ 7 x q(0), and 
the Q-ball charge depends linearly on the initial charge 
density. 

The second situation is the dependence on the angu- 
lar velocity (p'2{0) while <y5i(0) is fixed. We calculate for 
three different value of .^i(O): 10^ 10^, and 10^ In aU 
cases, linear dependence is still preserved when the ra- 
tio of </5i(0) and "^22(0) within two orders of magni- 
tude. However, if ^2(0) becomes smaller, the maximum 
Q-ball charge does not depend on the initial charge den- 
sity. This is due to the creation of the negative-charged 
Q-balls. The charge is determined only by iy9i(0). 

Negative charge Q-balls are formed when the (initial) 
angular velocity is rather small. Figure || shows an exam- 
ple. In this case, we see the largest Q-ball with positive 
charge, two large negative charge Q-balls, and one Q-ball 
with positive charge an order of magnitude smaller for 
four largest ones. Similar situations occur in the gauge 
mediation scenario , but the critical value of the ratio 
9?2(0)/(/?i(0) for the negative charge Q-ball formation is 
larger in the gravity mediation scenario. This is because 
the angular motion of the AD condensate is more circular 
and stable, and the produced Q-ball size is larger in the 
flatter potential, so that it is more difficult to reverse the 
angular velocity of the field within that size. 

In the actual situation, the AD field takes a very large 
vev before it rolls down to the origin of its potential, and 
the vev is determined by equating second and third terms 
in the potential (|]): 




FIG. 8. Configuration of positive and negative Q-balls 
on one-dimensional lattice. Here we take (pi{0) — 10^ and 
1/92(0) = lO'^. The four largest Q-balls have the charges (a) 
-2.9 X 10^^ (b) -5.6 X 10^^ (c) 6.8 x 10^^ and (d) 8.1 x 10^^ 
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The AD field begins to roll down when i? ~ m, so its am- 
plitude is (/J ~ (Ato)~^/^ ~ 2.4 x 10'' in the dimensionless 
parameters, where rh = m/M. At the same time, the AD 
field begins rotation because of the A-term of the form, 
VA-term ~ (\m / M)(f>'^ + h.c. If wc assumc that the initial 
angular velocity is the same order as the initial amplitude 
in the dimensionless units, we get the initial charge den- 
sity as q{0) = (pi(0)(/?^(0) - 6 X 10". We expect the lin- 
ear dependence between the initial charge density of the 
AD condensate and the produced largest Q-ball on three- 
dimensional lattices, as Qmax — x 10^. This is shown 
in Fig. ^, where we take such initial conditions as the lin- 
ear dependence is expected to hold, i.e., V'i(O) ~ Vii^)- 
Using this relation, we can estimate the maximum charge 
of the actually expected Q-balls is Qmax '^-^ 6 x 10^^. For 
the B-ball baryogenesis to work, the charge should exceed 
10^" !§]• Therefore, it may be a little difficult to reach 
this value in the parameters in the model. However, if 
we take X^cp^^/M^ instead of X'^cjj^/AP in the potential, 
as appears in the u'^(F(F flat direction the initial 

vacuum expectation value (vev) of the AD field is esti- 
mated as ip ^ (Ato'^)^^/'* ~ 7 X 10^". In this case, the 
initial AD charge density becomes ~ 5 x 10^^, and it im- 
plies that the maximum Q-ball charge reaches as large as 
~ 5 X 10^"^. Thus, we get enough value of the charge for 
B-ball baryogenesis. 



IV. B-BALL BARYOGENESIS AND ITS 
RESTRICTIONS TO THE PARTICLE PHYSICS 

As is known, baryon number and the amount of the 
dark matter can be directly related in the B-ball baryoge- 
nesis in the gravity- mediated SUSY breaking scenario [0 . 
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TABLE I. Fraction of the charge stored in Q-baUs for vari- 
ous values of K and (0)/<^i (0) on one-dimensional lattices. 
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FIG. 9. Dependence of charges on the initial charge den- 
sity g(0) = (^1(0)932(0) carried by the AD condensate on 
three-dimensional lattices. 



To this end, it is important to estimate how much charges 
are stored in the form of the Q-ball. In some cases, the 
fraction of the Q-ball charge may restrict the mass of the 
LSP, and vice versa |8|J|] . We have calculated for various 
initial conditions on one-, two-, and three-dimensional 
lattices, and find that almost all the charges are absorbed 
into Q-balls. This fact is also true when we take other 
values for parameters in the potential. In particular, we 
investigate for the fraction of Q-ball charge, changing K 
from —0.01 to —0.1. It was done by other method in 
Refs. [Q, and they concluded that when the absolute 
value I if I was larger, the less the fraction. However, our 
results differ from theirs. We collect them in Tables | and 
The former is the results from one-dimensional lat- 
tices with the box size iVA^ = 1024 x 0.02 = 20.48. The 
latter table shows the results calculated in three dimen- 
sions. In this case, the box size is iVA^ = 64 x 0.1 = 6.4. 
As can be seen, the fraction of the sum of charge of Q- 
balls to the total charge has no dependence on the value 
of K . Moreover, neither does it depend upon the ratio 
of </?i(0) and ^p^i^)- AH of them lead to a conclusion 
that almost all the charges are stored in Q-balls: that is. 

Following the argument of Refs. I^],!), the number den- 
sity of the baryon to that of the dark matter ratio can 
be written in terms of density parameters as 



rib 



b rUDM 



(22) 



where tun — ^ GeV is the nucleon mass. In the B-ball 
baryogenesis of the gravity-mediation scenario, B-balls 
decay into baryons and LSP neutralinos, so that the re- 
lation between the number density of baryon and dark 
matter is hum = NEfsn-b, where Nb is the number of 
neutralinos into which the AD field with a unit charge 
decays, and it is usually > 3. Here we assume no later 
annihilation of neutralinos. Using the conservative con- 
straint on the amount of the baryon number from the 
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99.1 


97.9 


98.6 


2.0 X 10'' 


99.0 


97.6 


98.3 


1.0 X 10" 


91.5 


97.5 


99.6 


8.0 X 10^ 


97.6 


95.5 


97.0 


6.0 X 10^ 


96.1 


97.4 


97.9 


4.0 X lO'^ 


99.4 


95.2 


99.7 


TABLE II. 


Fraction of the charg 


e stored in 


Q-balls for var- 


ious values of K and ip2(0)/v3i(0) 


on three-dimensional lat- 


tices. 








^2 \ K 


-0.01 


-0.05 


-0.1 


2.5 X lO'^ 


98.7% 


99.7% 


99.1% 


2.5 X 10" 


98.1 


99.4 


99.5 


2.5 X lO'^ 


98.4 


99.8 


99.2 



Big-Bang nucleosynthesis, 0.004 < flbh'^ < 0.023 111, we 
get a stringent constraint on the ncutralino mass 



7.1GeV<mv — ^ 



Nb 
3 



0.49 



Jb < 40.8GeV. 



(23) 



This bound is marginally consistent with fB~l and 
the accelerator experiment bounds such as > 24.2 
GeV [|l6|. Note that the constraint becomes more severe 
if Qdm is smaller than 1 as in the case, for example, 
that considerable fraction of the total energy density is 
stored in the form of the cosmological constant In 
this case, the annihilation of neutralinos must be taken 
place. 



V. MOVING Q-BALLS, THEIR INTERACTIONS, 
AND BREATHER-LIKE SOLITON 

As the consequence that the size of Q-balls is relatively 
small in the gravity-mediated SUSY breaking scenario, in 
a fixed volume, the coherent AD field breaks into larger 
numbers of Q-balls than in the gauge mediation scenario. 
Therefore, Q-balls can have somewhat large peculiar ve- 
locities, as opposed to Q-balls in gauge-mediation sce- 
nario. Actually, we observed moving Q-balls on the lat- 
tices in one, two, and three dimensions, but, unfortu- 
nately, Q-ball collisions (interactions) are observed only 
on one-dimensional lattices. This is not a surprise, since 
the impact parameter is small for small size Q-balls in two 
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or three dimensions. On the other hand, in one dimen- 
sion, Q-baUs must coUide if they have enough (initial) 
velocities. We see the following three patterns for the 
interactions: (a) passing through, (b) exchanging part of 
charges, and (c) merging. They are expressed symboli- 
cally as 

A + B — > B + A, (24a) 
A + B — >B' + A', (24b) 
A + B — yC. (24c) 

These situations are plotted in Fig. |o[ For the top three 
panels, they show the type (a), and two Q-balls with 
charges 4.0 x 10^^ and 1.8 x 10^^ are approaching, get to- 
gether with the charge 5.8 x 10^^, and finally pass through 
each other without changing their own charges. For the 
middle three panels, they represents the type (b). They 
exchange about 10% of their charges. In the bottom three 
panels, we show the merging process. 

Qualitatively, these processes can be divided by the 
relative velocity of two colliding Q-balls. If the relative 
velocity is large, they pass through each other without 
any (or negligible) charge exchange. When the velocity 
are smaller, two Q-balls exchange part of their charges. 
When the velocity is still slower, they merge into one, 
and it vibrates for a while. It can be a breather-like soli- 
ton, and an example is shown is Fig. ^ It repeats the 
double peaks and the single peak profiles just after the 
collision until it becomes stable state. During this pro- 
cess, we observed the decay of the charge by emitting 
very small Q-balls. For this particular example, about 
7% of its charge is lost until it finally becomes stable and 
conserves its charge from that time on. The decrease of 
charge can be explained also by the emission of scalar 
waves, but we cannot distinguish them in the resolution 
of our simulations. In addition to the merging process 
(c), we see a few inverse processes: the breaking into 
two. These three processes fa], (b), and (c) are very 
similar to the results of Ref . ||l8[ , where the collision of 
non-topological solitons for other type is studied numer- 
ically on two-dimensional lattices. Although we do not 
have a chance to see any collision in two or three dimen- 
sions, their properties may be very similar if it happens 
to occur. 



VI. CONCLUSIONS 

We have calculated the full non-linear dynamics of the 
complex scalar field, which represents some fiat direction 
carrying the baryonic charge in MSSM, in the context 
of the gravity-mediated SUSY breaking scenario. Since 
the scalar potential in this model is flatter than 0^, we 
have found that fluctuations develop and go non-linear to 
form non-topological solitons, Q-balls. As in the gauge- 
mediation scenario ||l^ , the size of a Q-ball is determined 
by the most amplifled mode, but this mode is completely 
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FIG. 10. Configurations of Q-balls for (a) passing through, 
(b) exchanging part of charges, and (c) merging. 
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FIG. 11. Configurations of merging Q-balls on 

one- dimensional lattices. Each of the panels show the time 
snapshots at from (1) r = 4.375 x 10** to (9) r = 4.775 x 10^ 
with the interval At = 0.05 x lO". 

determined by the model parameters m and K, and the 
size does not depend on the charge Q. On the other 
hand, the charge of Q-balls depends on the initial charge 
density of the Affleck-Dine field, and its dependence is 
linear. Therefore, large-charged Q-balls with relatively 
small size are formed in this scenario. 

Once Q-balls are formed, almost all the charges are 
absorbed into them in all the simulations we made, and 
only a tiny fraction of the charge is carried by the relic 
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AD field, but its amplitude is very small and fluctuates so 
that it may not be possible to regard it as a condensate. 
This leads to some interesting results. We can restrict the 
scenario of the baryogenesis, which has a direct relation 
to the amount of the dark matter, or the parameter in 
MSSM, such as the neutralino mass, can be constrained. 

We have also observed moving Q-balls, which is pecu- 
liar to the gravity-mediation scenario. In this case, larger 
numbers of Q-balls are formed in a fixed box size because 
of the relatively small Q-ball size, so the peculiar veloc- 
ities are larger than those in the gauge-mediation sce- 
nario. As a consequence, there are collisions of Q-balls. 
The probability of collision crucially depends on the spa- 
tial dimensionality, and we have not found any collision 
in two or three dimensions. We thus expect the proba- 
bility to be small in an actual situations. However, very 
interesting phenomena will occur, if collisions happen to 
take place. They are the charge exchange and merging 
to be large charge Q-balls. If the charge of a Q-ball be- 
comes larger, it will be more difficult to evaporate or to 
be dissociated. 
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